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The reduced matrix elements d2 and efe are computed in lattice QCD with Nf = 2 flavors of 
light dynamical (sea) quarks. For proton and neutron targets we obtain as our best estimates 



4 P) = 0.004(5) and d. 



(n) 



-0.001(3), respectively, in the MS scheme at Q — 5 GeV , while for 02 



we find — 0.077(12) and aS,' 1 ' = —0.005(5), where the errors are purely statistical. 
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I. INTRODUCTION 

The nucleon's second spin-dependent structure func- 
tion g 2 is of considerable phenomenological interest since 
at leading order in Q 2 it receives contributions from both 
twist- 2 and twist-3 operators. Consideration of 172 via the 
operator product expansion (OPE) offers the unique 
possibility of directly assessing higher-twist effects which 
go beyond a simple parton model interpretation. 

Neglecting quark masses and contributions of twist 
greater than two, one obtains the "Wandzura-Wilczek" 
relation Q 

g 2 (x, Q 2 ) » g¥ w (x, Q 2 ) = - 9l (x, Q 2 )+ f ^-g x {y, Q 2 ) , 

Jx y 

(i) 

depending only on the nucleon's first spin-dependent 
structure function, g\(x,Q 2 ). Including mass and gluon 
dependent terms up to and including twist-3, g 2 can be 
written Q 



From Eqs. QJ-Qi the moments of 172 are 



da; x n g 2 (x,Q 2 ) 



- dxx n 9l (x,Q 2 ) 
1 + 1 I Jo 



dx x 



n-l 



A I 



h T {x,Q 2 )+^x,Q 2 ) 



(4) 



A leading order OPE analysis with massless quarks 
shows that the moments of g\ and g 2 are given by Q] 



2 fdxx n gi {x,Q 2 ) = \ £ e[% 2 /Q 2 ,g(ri)a(JHv), 

(5) 



2 / dxx n g 2 (x,Q 2 ) 
Jo 



f—u,d 
1 n 
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g 2 (x,Q 2 )=g 2 vw {x,Q 2 )+ T2 {x,Q 2 ), 



(2) 



where 



92(x,Q 2 ) = - 



xd$P(jx) - e[% 2 IQ 2 ,g^))a^^)] , (6) 

for even n > for Eq. (JSJ and even n > 2 for Eq. JfjJ), 
where / runs over the light quark flavors and /1 denotes 
the renormalization scale. The reduced matrix elements 
a\P [fi) and d[P (p) are defined by Q 



'dj/A 



M 



h T (y,Q 2 ) + (i(y,Q 2 ) 



1 



-a 



if) 



.(3) 

The function }it(x,Q 2 ) denotes the transverse polariza- 
tion density and has twist two. The contribution from 
h T {x,Q 2 ) to <?2 is suppressed by the quark-to-nucleon 
mass ratio, m/M, and hence is small for physical up and 
down quarks. The twist-3 term £ arises from quark-gluon 
correlations. 
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• -D/^ ^ — traces . (9) 



Here D = D — D and e^, „ are the Wilson coefficients 
which depend on the ratio of scales n 2 /Q 2 , the running 
coupling constant and the quark charges 

e[% 2 /Q 2 ,g(ri) = Q (/)2 (i + 0(.9(m) 2 )) • (10) 

The symbol{- • • } ([•••]) indicates symmctrization (anti- 
symmetrization) of indices. The operator (JJJ has twist 
two, whereas the operator JBJ has twist three. Note that 
our definitions of a 2 and d 2 differ by a factor of two from 

those in HU. 

Using the equations of motion of massless QCD one 



-,5(/) 



such that 



can rewrite the twist- 3 operators O, r 

the dual gluon field strength tensor G^ v and the QCD 
coupling g appear. For n — 2 one finds 

°[l ( £dM = - f ^ (^mi7m2 + 7/n) V - traces, 

so we can define the reduced matrix element d 2 in the 
chiral limit also by (see, e.g., Ref. 0) 



traces 



traces \p, s) 

(12) 



This shows (setting fJ-i = fJ>2 = 0) that d 2 parametrizes 
the magnetic field component of the gluon field strength 
tensor which is parallel to the nucleon spin. Furthermore 
we have 



da; x£(x) 



(13) 



Hence, a calculation of d 2 (in the chiral limit) is especially 
interesting as it will provide insights into the size of the 
quark-gluon correlation term, £(x). 

The Wilson coefficients ijTU)) can be computed pertur- 

(f) ( f ) 

bativcly, while the reduced matrix elements a,/ and rf„ 
have to be computed non-perturbativcly. In the following 
we shall drop the flavor indices, unless they are necessary. 

A few years ago we computed the lowest non-trivial 
moment of g 2 in the quenched approximation Q. In 
this paper we give our results for the reduced matrix ele- 
ments a 2 and d 2 in full QCD, including Nj = 2 flavors of 
light dynamical (sea) quarks, using 0(a)-improved Wil- 
son fermions. We employ the same methods as in the 
quenched case, in particular the renormalization of the 
lattice operators is done entirely non-perturbatively. 



II. LATTICE OPERATORS AND 
RENORMALIZATION 

The lattice calculation divides into two separate tasks. 
The first task is to compute the nucleon matrix elements 



of the appropriate lattice operators. This was described 
in detail in The second task is to renormalize the 
operators. In the case of multiplicative renormalizability, 
the renormalized operator 0(/i) is related to the bare 
operator 0(a) by 



0(^)=Z (ap,)0(a), 



(14) 



where a is the lattice spacing. In our earlier work 
we computed the renormalization constants in pertur- 
bation theory to one-loop order. However, this does 
not account for mixing with lower-dimensional operators, 
which we encounter in the case of the reduced matrix 
elements dn ■ In Q an entirely non-perturbative solu- 
tion to this problem was presented for quenched lattice 
QCD. Here we shall apply the same approach. We im- 
pose the (MOM-like) renormalization condition [ToL ITU 
(which can also be used in the continuum) 

±Tr( q (p)\0(n)\ q (p))\( q (p)\0(a)\ q (p))r 

L j p'=fi 

(15) 

where \q{p)} is a quark state of momentum p in Landau 
gauge. 

In the following we shall restrict ourselves to the case 
n = 2. Furthermore, we consider quark-line connected 
diagrams only, as calculations of quark-line disconnected 
diagrams are extremely computationally expensive. In an 
attempt to improve on our earlier analysis 0, we simu- 
late with two non-vanishing values for the nucleon mo- 
mentum, pi = (p, 0,0) and p 2 = (0,p, 0), together with 
two different polarization directions, described by the 
matrices Y x = \(\ + 74) i7 5 7i and T 2 = |(1 + 74) i7572- 
Here p = 2tt/Ls denotes the smallest non-zero momen- 
tum available on a periodic lattice of spatial extent Ls- 
We consider the two combinations px/T^ and P2/T1. For 
the twist-2 matrix element a 2 we use in both cases the 
operator 



)5 _. q{5} 
{214} ' ^ 



(16) 
as in ■ 

For the twist-3 matrix element d 2 we need to use differ- 
ent operators for our two momentum/polarization com- 
binations. Forpi/r2 andp2/ri we take 
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respectively. In the following we shall suppress the index 
of 0^ unless it is needed. The operators and 
belong to the representations 7^ and , respectively, 
of the hypercubic group H(4) [Tj]. The operator has 
dimension five and C-parity +. It turns out that there 
exist two operators of dimension four and five, respec- 
tively, transforming identically under H{4) and having 
the same C-parity, with which can mix: 



<7i 3 £>i - a 43 D 4 i> =■ O", (19) 



- l4 D 3 D 4 + l4 D 4 D^jij=:O , (20) 

for P1/T2, and similarly for P2/T1 with 1 — > 2. We use 
the definition er M „ = (i/2)^, 

The operator (|20() mixes with with a coefficient of 
order g 2 and vanishes in the tree approximation between 
quark states. We therefore neglect its contribution to the 
renormalization of . The operator O a , on the other 
hand, contributes with a coefficient cx a~ x and hence 
must be kept. We then remain with 



= Z^{a^)0 [5] {a) + -Z a ( af i)O a (a). 



(21) 



The renormalization constant and the mixing coeffi- 
cient Z a are determined from 



Tr (q(p)\OW(ri\q(p)) [(q(p)\0 [5] (o)|?(p)) | trc 
Tr(q(jp)\OW(n)\q(p)) \(q(p)\0 ° (a)\q(p)) | trc 



(22) 

-1 

= 0. 

p 2 —fj, 2 

(23) 



Rewriting Eq. 1(21)1 as 

0^) = Z^){0^ia) + \^^0%a)^ , (24) 

we see that O^(fx) will have a multiplicative dependence 
on fi only if the ratio Z° '(afj,)/Z^(an) does not depend 
on /i, which should happen for large enough values of the 
renormalization scale. The scale dependence will then 
completely reside in Z^ . 



III. SIMULATION DETAILS 

To reduce cut-off effects, we use non-perturbatively 
0(a) improved Wilson fermions. The calculation is done 
at four different values of the coupling, /?, and at three 
different sea quark masses each. The latter are specified 
by the hopping parameter K soa . We use the force param- 
eter ro to set the scale, with tq — 0.467 feu. Our lattice 
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FIG. 1: Parameters of our dynamical gauge field configura- 
tions, together with lines of constant ro/a (solid lines) and 
lines of constant mps ro (dashed lines). The simulations are 
done on 24 3 x 48 (x) and 16 3 x 32 (O) lattices, respectively. 



TABLE I: Lattice parameters: gauge coupling /3, sea quark 
hopping parameter K S ca, lattice volume, number of trajecto- 



ries, 


ro/a and pseudoscalar meson 


mass. 







^sca 


Volume 




ro/a 


mpso 


5.20 


0.13420 


16 3 


x 32 


0(5000) 


4.077(70) 


0.5847(12) 


5.20 


0.13500 


16 3 


x 32 


O(8000) 


4.754(45) 


0.4148(13) 


5.20 


0.13550 


16 3 


x 32 


O(8000) 


5.041(53) 


0.2907(15) 


5.25 


0.13460 


16 3 


x 32 


0(5800) 


4.737(50) 


0.4932(10) 


5.25 


0.13520 


16 3 


x 32 


O(8000) 


5.138(55) 


0.3821(13) 


5.25 


0.13575 


24 3 


x 48 


0(5900) 


5.532(40) 


0.25638(70) 


5.29 


0.13400 


16 3 


x 32 


O(4000) 


4.813(82) 


0.5767(11) 


5.29 


0.13500 


16 3 


x 32 


0(5600) 


5.227(75) 


0.42057(92) 


5.29 


0.13550 


24 3 


x 48 


O(2000) 


5.566(64) 


0.32688(70) 


5.40 


0.13500 


24 3 


x 48 


0(3700) 


6.092(67) 


0.40301(43) 


5.40 


0.13560 


24 3 


x 48 


0(3500) 


6.381(53) 


0.31232(67) 


5.40 


0.13610 


24 3 


x 48 


0(3500) 


6.714(64) 


0.22120(80) 



spacings range from a = 0.07 to 0.09 fm. The actual 
parameters, as well as the corresponding values of ro/a 
and the pseudoscalar meson masses, are given in Table H] 
and shown pictorially in Fig. ^ 

The quark matrix elements for the renormalization 
constants are computed using a momentum source |l II . 
Performing the Fourier transform at the source sup- 
presses the effect of fluctuations: The statistical error in 
this case is cx (T^^ CO nf) _1 ^ 2 for A^onf configurations on 
a lattice of volume V, resulting in small statistical un- 
certainties even for a small number of configurations, at 
least five in our case. Hence, the main source of statisti- 
cal uncertainty in our final results is from the calculation 
of the bare matrix elements, not the Z values. 

Nuclcon matrix elements are determined from the ratio 
of three-point to two-point correlation functions 



n{t,r-p-o) 



Cr(t,T;p,Q) 
C 2 (t,p) 



(25) 



4 



where C2 is the vmpolarized baryon two-point function 
with a source at time and sink at time t, while the three- 
point function Cr has an operator O insertion at time r. 
To improve our signal for non-zero momentum we average 
over both polarization/momentum combinations. 

Correlation functions are calculated on configurations 
taken at a distance of 5-10 trajectories using 4-8 different 
locations of the fermion source. We use binning to obtain 
an effective distance of 20 trajectories. The size of the 
bins has little effect on the error, which indicates auto- 
correlations are small. 



IV. COMPUTATION OF RENORMALIZATION 
CONSTANTS 




10 20 30 40 50 60 70 
ff (GeV 2 ) 



The twist-2 operator defined in Eq. (|16[) is renor- 
malized multiplicatively with the renormalization factor 
Z^(afi), while the renormalization of the twist-3 oper- 
ators in Eqs. 1|17|) . (|18|l is more complicated due to the 
mixing effects described in Section [n] Since the renor- 
malization of o[ 5 ' and is identical (up to lattice arte- 
facts) we consider only of\ 

The calculation of the non-perturbative renormaliza- 
tion factors is a non-trivial exercise, the full details of 
which are beyond the scope of this paper. Here we re- 
strict ourselves to a short outline of the procedure. More 
details can be found in Section 5.2.3 of Ref. and a 
fuller account will be given in a forthcoming publication. 

Firstly, a chiral extrapolation of the non-perturbative 
renormalization factors is performed at fixed (3 and fixed 
momentum. The extrapolation is performed linearly in 
(romps) 2 = ((r / a)amps) 2 , where for each value of (3 we 
use the chirally extrapolated value of r /a (see Table 3 
of Ref. ^3)- We then apply continuum perturbation 
theory to calculate the renormalization group invariant 
renormalization factor Zrgi from the chirally extrapo- 
lated Zs |13J. This can be done in various schemes, 
e.g., the MS scheme, and should lead for any scheme 
to the same momentum- independent value of Zrgi, at 
least for sufficiently large momenta. For this step, we use 
r o^MS = 0-617 In Fig. [3 we show the /Lt-dependence 
of .Zj^qj computed in the MS scheme and in a continuum 
MOM scheme at (3 = 5.40. While in both cases a reason- 
able plateau appears, the plateau values do not coincide 
exactly, and we take the difference as a measure of the 
uncertainty of our Zs, caused by our incomplete knowl- 
edge of the perturbative expansion. 

The final step requires -Zrgi to be converted to Zy^ at 
some renormalization scale, which is done perturbatively, 
and the result depends on the value of A^jg in physical 
units. From rpAj^g = 0.617 and r$ = 0.467 fm we obtain 



A- 



MS 



261 MeV. 



As mentioned above, the renormalization of the twist- 
3 operator in Eqs. I|17|) . I|18|) has further complications 
due to the mixing effects described in Section [D] In this 
case it is unclear how to convert our MOM results to the 



FIG. 2: -Z^cji calculated in the MS scheme (circles) and in a 
MOM scheme (filled squares) at j3 = 5.40. The scale is fixed 
using ro = 0.467fm. 
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FIG. 3: The ratio Z" 7 (afi)/Z [5] (afj.) at (3 = 5.40 



MS scheme. So we shall stick to the MOM numbers. For 
the comparison of our results with experimental determi- 
nations this does not cause problems, because no QCD 
corrections have been taken into account in the analysis 
of the experiments and hence different schemes are not 
distinguished. 

In Fig. |2| we plot the ratio Z a (a/j,)/Z^ (afj.) as a func- 
tion of /i for (3 — 5.40. As expected, a plateau develops 
for larger values of /i, and therefore the operator (/i) 
only depends on \x multiplicatively. 
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RESULTS FOR REDUCED MATRIX 
ELEMENTS 



0.15 



In order to compute the reduced matrix elements in 
Eqs. and (JSJ, we calculate the ratio of three- to two- 
point correlation functions 1Z, as given in Eq. Q25[l. for the 
operators defined in Eqs. (|16[1 -H20 [1 . The bare operator 
matrix elements are obtained from the ratio 1Z by 



-Mpa 2 , TZd 2 



— \hl P d 2 . (26) 



We define the continuum quark fields by ^/2k sc& times 
the lattice quark fields. The factor for 72.^2 is the same 
for all three operators , O a and O . 

In Tables ITT1 and Hill we present our results for the bare 
matrix elements of the operators 0^ 5 \ O^, O a and O 
defined in Eqs. ^SJ-|J2DJ for u and d quarks in the proton. 

The corresponding renormalized (reduced) matrix el- 
ements for the renormalization scale /i 2 = 5 GeV 2 are 
given in Tables IIVI and [V] While the superscripts (it) 
and (d) again refer to u and d quarks in the proton, the 
matrix elements for proton and neutron targets are de- 
noted by (p) and (n), respectively For a 2 the latter are 
given by 



-0) 
h 

» 



)(u)2 q («) 
,(«) 



= Q(<0 2 , 



A d ) 



q(«) 2 , 



and similarly for d 2 . The renormalized values of d 2 
f = u, d in the proton are calculated from 



(/) 



(27) 
(28) 

for 
(29) 



In the lines for K sca = k c , Tables IIVI and |V| contain 
results in the chiral limit, obtained by an extrapola- 
tion linear in (romps) 2 - The scale has been fixed from 
the value of ro/a at the respective quark masses using 
ro = 0.467 fm. Alternatively, we could have worked with 
the chirally extrapolated values of ro/a. This would in- 



4 /) = ^ 51 4 5](/) + -^rf? (/) 
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FIG. 4: The chirally extrapolated reduced matrix element a 2 
for the proton target renormalized at the scale fj, 2 = Q 2 — 5 
GeV 2 as a function of the lattice spacing a. The crosses denote 
phenomenological determinations. 
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crease d 2 p) and d> 2 > by up to twice the statistical error FIG. 5: The chirally extrapolated reduced matrix element a 2 
but would leave the other observables almost unaffected. for the neutron target renormalized at the scale fi 2 = Q 2 = 5 
On the other hand, setting ro = 0.5 fm or varying ro A^g- GeV 2 as a function of the lattice spacing a. The cross denotes 
between 0.572 and 0.662 (corresponding to the combined the phenomenological value, 
statistical and systematic errors given in Ref. Q) leads 
only to rather small changes in the final results. 

Let us first focus on the results for the twist-2 ma- 
trix element a 2 ■ In Fig. 0] we show the chirally extrap- 
olated renormalized results for a 2 in the proton in the 
MS scheme as a function of the lattice spacing a. It 
should however be noted that the data at (3 = 5.20, i.e., 
those for the largest lattice spacing are to be considered 
with caution, because potentially they are affected by lat- 
tice artefacts. For a 2 the dependence on the quark mass 
turns out to be rather small. On the other hand, we do 
not attempt a continuum extrapolation of the chirally 
extrapolated results. Instead we take the value at our 
smallest lattice spacing ((3 = 5.4) as our best estimate: 



Au) 



= 0.077(12). This is consistent with earlier quenched 
results |7J, indicating that quenching effects are small. 

At the physical pion mass, we compare with two re- 
sults taken from the literature which are obtained from 
an analysis of experimental data. The larger value is 
taken from an earlier analysis performed by Abe et al. 
[4j , while the lower point is extracted from a recent anal- 
ysis by Osipenko et al. 0] with the help of the per- 
turbative Wilson coefficient. In the MS scheme with an- 
ticommuting 75, we use the two-loo p ex pression for the 
Wilson coefficient described in Ref. |20|. To avoid large 
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TABLE II: Bare (unrenormalized) matrix elements 02, d 2 , for u and d quarks in the proton for our entire set of ((3, K SC a) 
combinations. 



p 


^sca 


4 


Id) 
4 


4 


j[5] (d) 


5.20 


0.13420 


0.142(18) 


-0.0318(78) 


-0.0143(23) 


0.0005(14) 


5.20 


0.13500 


0.123(22) 


-0.032(11) 


-0.0329(59) 


0.0094(35) 


5.20 


0.13550 


0.131(32) 


-0.061(22) 


-0.057(14) 


0.0064(59) 


5.25 


0.13460 


0.113(12) 


-0.0389(51) 


-0.0165(25) 


0.0023(13) 


5.25 


0.13520 


0.110(19) 


-0.0281(74) 


-0.0310(39) 


0.0069(17) 


5.25 


0.13575 


0.1107(74) 


-0.0345(47) 


-0.0575(28) 


0.0074(15) 


5.29 


0.13400 


0.1141(77) 


-0.0255(35) 


-0.0033(11) 


-0.00009(63) 


5.29 


0.13500 


0.0989(90) 


-0.0281(45) 


-0.0252(19) 


0.0046(11) 


5.29 


0.13550 


0.1228(65) 


-0.0302(26) 


-0.0468(23) 


0.00783(92) 


5.40 


0.13500 


0.1195(44) 


-0.0227(24) 


-0.02135(99) 


0.00232(61) 


5.40 


0.13560 


0.1238(63) 


-0.0331(34) 


-0.0445(26) 


0.0069(11) 


5.40 


0.13610 


0.127(13) 


-0.0277(60) 


-0.0674(48) 


0.0103(25) 



TABLE III: Bare (unrenormalized) matrix elements d%/a and d® for u and d quarks in the proton for our entire set of (/3, ft sca ) 
combinations. 



p 




d 2 /a 


dl W /a 


4 M 




5.20 


0.13420 


-0.220(19) 


0.046(8) 


-0.0312(46) 


0.0096(22) 


5.20 


0.13500 


-0.305(29) 


0.077(13) 


-0.039(10) 


0.0145(49) 


5.20 


0.13550 


-0.395(60) 


0.080(21) 


-0.063(14) 


0.0194(75) 


5.25 


0.13460 


-0.252(17) 


0.045(6) 


-0.0371(34) 


0.0150(28) 


5.25 


0.13520 


-0.239(23) 


0.063(10) 


-0.0329(61) 


0.0131(42) 


5.25 


0.13575 


-0.353(13) 


0.0638(44) 


-0.0463(39) 


0.0141(20) 


5.29 


0.13400 


-0.213(9) 


0.0379(35) 


-0.0322(23) 


0.0086(12) 


5.29 


0.13500 


-0.258(13) 


0.0518(42) 


-0.0312(34) 


0.0118(21) 


5.29 


0.13550 


-0.338(10) 


0.0651(36) 


-0.0390(25) 


0.0120(13) 


5.40 


0.13500 


-0.301(8) 


0.0595(33) 


-0.0396(18) 


0.01231(84) 


5.40 


0.13560 


-0.385(15) 


0.0723(50) 


-0.0502(26) 


0.0137(15) 


5.40 


0.13610 


-0.420(25) 


0.087(9) 


-0.0411(60) 


0.0178(39) 



logarithms, we set Q 2 — fi 2 — 5 GeV 2 to obtain 



proposed in 



,(/) _ (f)2 



x 1.03075. 



(30) 



We do not see exact agreement between our chirally ex- 
trapolated value and those obtained from experimental 
data, but there are still several sources of systematic error 
in our final number. Firstly, our simulation only involves 
the calculation of connected quark diagrams. That is, 
we do not consider the (computationally expensive) case 
where an operator couples to a disconnected quark loop, 
although such disconnected diagrams are not expected 
to contribute in the large x region. Secondly, our results 
are restricted to the heavy pion world, mps > 550 MeV. 
In this region we observe a linear dependence of our re- 
sults on TOpg. A more advanced functional form guided 
by chiral perturbation theory, such as those proposed for 
the moments of unpolarized nucleon structure functions 
or nucleon magnetic moments ^7]]> ma y be required. 
One such form has been suggested in 0, but only for 
iso- vector matrix elements. So we attempt to gain an 
estimate of the systematic uncertainty due to our linear 

extrapolation by comparing results for a 2 u ^ in the chi- 
ral limit using both a linear extrapolation and the form 



(u-d) 



a. 1 



a, I 1 + C L NA m; 1o S -^TT 



+b 2 



(31) 



where the authors recommend a preferred value for the 
LNA coefficient as c L na = -(0.48^ + l)/(47r/ 7r ) 2 and b 2 
is constrained by the heavy quark limit to be 



Au-o 



5 

27 



(u-d) 



(1 - ^ 2 c LNA ) 



(32) 



We set 11 = 0.25 GeV as proposed in 0] and find at (3 = 
0.214(29) employing a linear extrapolation 
0.183(9) using Eq. ffiTfr. suggesting there is 



5.29, a 



(u-d) _ 



(u-d) 



and a. 

a 15% systematic error in our linear extrapolation. 

Finally, we have not considered finite size effects 19] 
in this work, and our data do not yet allow us to perform 
a decent continuum extrapolation. 

Our results for a 2 in the neutron are shown in Fig. [S] 
They are hardly different from zero. Taking again the 
value for (3 = 5.4 as our best estimate, we end up with 
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TABLE IV: Renormalized matrix elements for the renormalization scale /i 2 = 5 GeV 2 in the MS scheme. The superscripts (u) 
and (d) refer to u and d quarks in the proton. 



p 


^sca 


(u) 






4 d) 


5.20 


0.13420 


0.194(27) 


-0.044(11) 


0.0360(59) 


-0.0113(29) 


5.20 


0.13500 


0.168(32) 


-0.044(15) 


0.039(12) 


-0.0082(65) 


5.20 


0.13550 


0.179(45) 


-0.083(30) 


0.034(28) 


-0.015(11) 


5.20 




0.154(65) 


-0.079(37) 


0.040(31) 


-0.011(14) 


5.25 


0.13460 


0.154(19) 


-0.0532(76) 


0.0335(53) 


-0.0070(24) 


5.25 


0.13520 


0.150(27) 


-0.038(10) 


0.0109(79) 


-0.0047(34) 


5.25 


0.13575 


0.151(13) 


-0.0472(70) 


0.0024(54) 


-0.0050(25) 


5.25 


Kc 


0.149(24) 


-0.042(12) 


-0.0169(89) 


-0.0036(41) 


5.29 


0.13400 


0.159(14) 


-0.0356(53) 


0.0468(27) 


-0.0094(13) 


5.29 


0.13500 


0.138(15) 


-0.0392(67) 


0.0284(43) 


-0.0064(20) 


5.29 


0.13550 


0.171(13) 


-0.0421(43) 


0.0201(44) 


-0.0056(17) 


5.29 




0.167(24) 


-0.0469(84) 


-0.0008(70) 


-0.0026(28) 


5.40 


0.13500 


0.170(12) 


-0.0323(39) 


0.0499(27) 


-0.0127(13) 


5.40 


0.13560 


0.176(13) 


-0.0471(55) 


0.0401(57) 


-0.0097(22) 


5.40 


0.13610 


0.181(21) 


-0.0394(88) 


0.019(10) 


-0.0094(46) 


5.40 


/> c 


0.187(28) 


-0.056(11) 


0.010(12) 


-0.0056(50) 



TABLE V: Renormalized matrix elements for the renormalization scale fj, 2 = 5 GeV 2 in the MS scheme. The superscripts (p) 
and (n) denote the matrix elements for proton and neutron targets, respectively. 



p 


^sca 


4 P) 


(n) 


4 P) 


4 n) 


5.20 


0.13420 


0.081(12) 


0.0022(55) 


0.0148(26) 


-0.0010(14) 


5.20 


0.13500 


0.070(14) 


-0.0008(75) 


0.0166(55) 


0.0008(32) 


5.20 


0.13550 


0.070(20) 


-0.017(14) 


0.013(13) 


-0.0028(58) 


5.20 


/if 


0.058(29) 


-0.020(18) 


0.017(14) 


-0.0002(71) 


5.25 


0.13460 


0.0627(82) 


-0.0065(36) 


0.0141(24) 


0.0006(12) 


5.25 


0.13520 


0.063(12) 


-0.0004(53) 


0.0043(36) 


-0.0009(18) 


5.25 


0.13575 


0.0620(58) 


-0.0041(31) 


0.0005(24) 


-0.0019(13) 


5.25 


t\ c 


0.062(10) 


-0.0024(53) 


-0.0079(40) 


-0.0035(21) 


5.29 


0.13400 


0.0668(61) 


0.0019(25) 


0.0198(12) 


0.00105(64) 


5.29 


0.13500 


0.0570(65) 


-0.0021(31) 


0.0119(19) 


0.00031(99) 


5.29 


0.13550 


0.0715(57) 


0.0003(19) 


0.0083(20) 


-0.00028(89) 


5.29 


ft" C 


0.069(10) 


-0.0015(38) 


-0.0006(31) 


-0.0012(15) 


5.40 


0.13500 


0.0720(50) 


0.0045(17) 


0.0208(12) 


-0.00009(63) 


5.40 


0.13560 


0.0731(58) 


-0.0014(24) 


0.0168(25) 


0.0001(11) 


5.40 


0.13610 


0.0760(93) 


0.0026(43) 


0.0072(46) 


-0.0021(23) 


5.40 




0.077(12) 


-0.0048(53) 


0.0039(54) 


-0.0013(26) 



a 2 = —0.005(5), in agreement with the result from the 
analysis of Abe et al. |j| . 

From and in the chiral limit we calculate (see 
Eq. (JjJJ) the second moment of the polarized structure 
function g± for the proton and neutron. Using the Wilson 
coefficient given in Eq. 1|30[) we find 

/ dxx 2 g p 1 (x,Q 2 ) = — a p 2 = 0.0170(18) , (33) 
Jo 4 

/ dxx 2 g^(x,Q 2 ) = im ° 75 a P = -0.0013(8). (34) 
Jo 4 

We now turn our attention to the second moment of 
52- We find that our data for <i 2 also exhibit a linear 
behavior in mp S . While this is not unexpected at the 
large pion masses where our simulations are performed, 



this linear behavior will not necessarily continue near the 
chiral limit. Unfortunately, the dependence of on the 
pion mass near the chiral limit is not yet known. There- 
fore in this work we perform only a linear extrapolation 
of c?2 to the chiral limit. In Figs. and [7| we plot some 
of the data versus (romps) 2 together with the linear ex- 
trapolations. The chirally extrapolated results for d,2 in 
the proton and neutron are shown in Figs. |H] and re- 
spectively. At our smallest lattice spacing we obtain in 
the chiral limit 

4 P) = 0.004(5), (35) 
4 n) = -0.001(3). (36) 

The errors are statistical only. Taking the behavior of 
&2 as a guide, the chiral extrapolation might intro- 
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i 

0.008 



0.01 



FIG. 6: The chiral extrapolation of the reduced matrix el- 
ement d,2 for the proton target renormalized at the scale 



5 GeV^ 



FIG. 8: The chirally extrapolated reduced matrix element d,2 
for the proton target renormalized at the scale fj, 2 = Q 2 — 5 
GeV 2 as a function of the lattice spacing a. The cross denotes 
the phenomenological value. 



0.008 



0.004 - 




-0.008 - 



0.04 
0.03 
0.02 



= f 0.01 



0.0 - 
-0.01 - 
-0.02 



(r mJ 



0.0 0.002 0.004 0.006 0.008 0.01 
a 2 (fm 2 ) 



FIG. 7: The chiral extrapolation of the reduced matrix el- 
ement d,2 for the neutron target renormalized at the scale 



/'" 



= 5 GeV 



FIG. 9: The chirally extrapolated reduced matrix element d,2 
for the neutron target renormalized at the scale fi 2 = Q 2 — 5 
GeV 2 as a function of the lattice spacing a. The cross denotes 
the phenomenological value. 



u 2 = 

6 



/ dxx 2 q 2 (x) = f dxx- 
Jo K ' Jo 3 



ducc a 15% systematic uncertainty. For the other 

systematic uncertainties discussed above would amount explicit Q 2 dependence) 

(n) 

to an additional error of about 0.005, while d 2 is almost 
unaffected. Our result for the proton agrees very well 
with the experimental number |5j, while for the neutron 
the experimental result differs from ours by two standard 
deviations. A more precise experimental value would be 
most desirable in case of the neutron. 

From Eq. J1J, the moments of g 2 receive contributions 
from gi and ^2, the second of which contains a mass 
dependent term and a gluon insertion dependent term. 
From Eq. (jSJ) , the second moment of 52 is (dropping the 



— hr(x)+t(x) 



(37) 

so if d 2 vanishes in the chiral limit, then J Q dxx£(x) must 
also vanish. Our results lead us to conclude that for the 
n = 2 moment the Wandzura-Wilczek relation 



/ dxx 2 g 2 (x,Q 2 ) = ~ I dxx 2 g!(x,Q' 
Jo a Jo 



) (38) 
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is satisfied within errors for both proton and neutron tar- 
gets. 

From the expression in Eq. (JSJ), we also expect the first 
moment of to vanish in the chiral limit. Combining 
these two observations with the Burkhardt-Cottingham 
sum rule |2l|, J* Q2(x)dx = 0, and the knowledge that 
from elastic scattering processes gi receives non-trivial 
higher-twist contributions at x = 1 (see, for example, 
Eqs. (4), (5) of H3)i we expect that there should be 
some sort of smooth transition at intermediate x, which 
presents an interesting challenge for the planned experi- 
ments at JLab 221. 



VI. CONCLUSIONS 

We have calculated the second moments of the proton 
and neutron's spin-dependent g\ and g2 structure func- 
tions in lattice QCD with two flavors of C(a)-improved 
Wilson fcrmions. A key feature of our investigation is 
the use of non-perturbative renormalization and the in- 
clusion of operator mixing in our extraction of the twist-2 
and twist-3 matrix elements. 

Our result for a 2 p ^ = 0.077(12) for the proton is some- 
what larger than what follows from analyses of experi- 
mental data, while for the corresponding result for the 



neutron, we find a small but negative value, a 2 = 
—0.005(5), in agreement with experiment. Note that the 
errors are purely statistical and do not include any sys- 
tematic uncertainties, although we estimate a systematic 
uncertainty of approximately 15% arising from the chiral 
extrapolation. 

For the twist-3 matrix element, d,2, our results agree 
very well with experiment and are consistent with zero, 
leading us to the conclusion that higher-twist effects oc- 
cur only at large or intermediate x. 
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